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We present a systematic study of the dijet suppression at RHIC using the VNI/BMS parton 
cascade. We examine the modification of the dijet asymmetry Aj and the within-cone transverse 
energy distribution (jet-shape) along with partonic fragmentation distributions z and jt in terms of: 
q; the path length of leading and sub-leading jets; cuts on the jet energy distributions; jet cone angle 
and the jet-medium interaction mechanism. We find that Aj is most sensitive to q and relatively 
insensitive to the nature of the jet-medium interaction mechanism. The jet profile is dominated by 
q and the nature of the interaction mechanism. The partonic fragmentation distributions clearly 
show the jet modification and differentiate between elastic and radiative+elastic modes. 



I. INTRODUCTION 

Recent measurements showing the strong modification 
of high energy E t ~ 100 - 200 GeV dijets at the LHC 
[TJ [2] in Pb+Pb collisions have sparked an interest in us- 
ing dijets to quantify jet modification in hot QCD mat- 
ter. Proposed upgrades to the PHENIX experiment, and 
continuing development of jet reconstruction at STAR, 
promise access to dijets in a different kinematic region, 
E t ~ 15 - 65 GeV at RHIC. The dijet asymmetry 



where E tt e is the transverse energy of the leading jet and 
E t s is that of the sub-leading jet, has been successfully 
reproduced by various authors with a variety of models 
[3HB]. Recent work by Renk [7] has indicated that these 
results may be relatively insensitive to the fine details of 
jet energy loss, suggesting that Aj may be unsuitable for 
tomographic purposes. If this is indeed the case are there 
other more differential dijet observables for jet measure- 
ments at RHIC? 

With these issues in mind we have used the VNI/BMS 
parton cascade to systematically explore dijet suppres- 
sion at RHIC energy scales under controlled medium 
conditions. We examine the modification of dijets un- 
der variation of the medium radius and temperature, the 
strong coupling constant and the jet definition in terms 
of the Anti-Kt cone angle and the leading jet energy cut. 

The VNI/BMS parton cascade model [H [H] provides 
access to the full jet/medium development at a fixed q. 
We run the code in a static uniform-medium mode. Here 
the medium is modeled as torus of a given radius. The di- 
jet propagation lengths are generated as chords centered 
on uniformly sampled hard-collision vertices within this 
torus. The model includes a partonic medium which is 
treated on an equal footing with the jet. This allows jet 
partons to escape into the medium and vice-versa. In 
keeping with the infinite static model of the medium we 
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do not consider hadronization of the jet. As such all re- 
sults are presented at the partonic level only. Finally, we 
use anti-kt jet reconstruction throughout to give a some- 
what realistic treatment of the jet measurement process. 

We begin with a brief description of our model. We 
then show results for the leading parton energy loss and 
discuss the modification of dijets at RHIC in terms of: 
the asymmetry Aj, the jet radial profile, and partonic 
fragmentation distributions z and j t . 

II. THE PARTON CASCADE MODEL 

The parton cascade model (PCM) is a Monte-Carlo 
implementation of the relativistic Boltzmann transport 
of quarks and gluons 

P ^F k (x,p) =J2CiF k (x,p). (2) 

i 

The collision term Ct includes all possible 2 — !> 2 interac- 
tions and final-state radiation 1 — > n 

C t F k (x,p) = ^ f HdTjlMifx 

1 j 

S 4 (P in -P out )D(F k (x,p}), (3) 

dTj is the Lorentz invariant phase space for the process j, 
D is the collision flux factor and Si is a process dependent 
normalization factor. A geometric interpretation of the 
total cross-section is used to select pairs of partons for 
interaction. Between collisions, the partons propagate 
along straight line trajectories. 

In the VNI/BMS PCM outgoing off-shell partons are 
brought back on-shell through a medium modified time- 
like branching. The partons created in this process is 
subject to a Monte-Carlo LPM effect, see below. 

The strong coupling constant for scatterings is nomi- 
nally held fixed at a s = 0.3 although we shall explore the 
effects of its variation below. A QGP medium is simu- 
lated as a box of thermal quarks and gluons generated at 
some fixed temperature. Periodic boundary conditions 
are imposed on the box, whose size is selected to be large 
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enough that if a simulated jet wraps around it will not 
interact with its own tail. The jet transport coefficient 
characterizing the amount of transverse momentum ac- 
quired by particles in transit through a perturbative ther- 
mal medium q can be analytically computed [101 1 1 1 j and 
also deduced from the simulation. The analytical result 
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where T is the medium temperature mo (T) is the Debye 
screening mass, J\f is the medium density and q ma x is 
a cutoff resulting from the parton scattering kinematics. 
See Fig.[l]for a comparison of results from VNI/BMS and 
Q. The partonic contents of jets created by a suitable 
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FIG. 1: The transport coefficient q for quark jets shown 
as a function of jet energy and medium temperature, 
the VNI/BMS data (crosses and open-triangles) agrees 
well with the theoretical expression (solid lines) Q. 

event generator and evolved down to Qq — 1 GeV, in this 
study PYTHIA 8 QH EES], are injected into the medium 
one at a time, they propagate for a given distance, and 
their evolution is recorded. The jet evolution is recorded 
at run-time without reference to any particular jet def- 
inition (cone-angle, energy cuts, jet reconstruction algo- 
rithm, etc). Instead each particle in the initially inserted 
jet is marked as being "jetty" , this jettiness tag is then 
iteratively applied to all partons that interact with an 
already jetty-labeled parton. The entire time-evolution 
of the jet can then be reconstructed using a suitable jet 
finder (Fast Jet [14]) and jet definition. In this way the 
influence of varying jet definitions can be readily stud- 
ied. We emphasize that the medium partons propagate 
along with the jet during this process, continually inter- 



acting and mixing with the evolving jet. This medium 
back-reaction is a unique feature of parton cascades. 

VNI/BMS is a simple enough jet-suppression model. 
While it may not be suitable for use as a full event gener- 
ator, lacking hydrodynamical flow and hadronization, the 
active medium and medium-modified radiation make it a 
useful test bed for many but not all jet modification fea- 
tures. As we shall show below the effects of hadronization 
are rather important for certain intra-jet observables. We 
shall reexamine these observables in a future work with 
a hadronization solution implemented in the model. 



III. A MONTE-CARLO LPM EFFECT 

The emission of radiation in QCD is not an instan- 
taneous process. There is a time period during which 
the radiated quanta and the radiator are in a correlated 
state. If this takes place in a dense medium the corre- 
lated wave function may interact with the medium con- 
stituents leading to a modification of the decoherence 
process. If the emitting parton interacts with several 
additional medium scattering centers then quantum in- 
terference effects lead to an in medium path-length de- 
pendence for the amount of energy radiated, this is the 
Landau Pomeranchuk Migdal (LPM) effect [T51[T5] . 

The formation time of a parton with energy E and 
time-like virtuality Q 
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is the lifetime of this virtual correlated state and charac- 
terizes the emission process. 

We implement an approximation to the full pQCD de- 
scription of the interaction of the radiating system with 
the medium using a local Monte-Carlo routine in the 
style of Zapp and Wiedemann [17] . This prescription 
reproduces the leading BDMPS-Z [HI US] result for light- 
parton energy loss in a QGP medium AE ~ 1? . This 
method is particularly appealing since it requires no ar- 
tificial parametrization of the radiative process, it is a 
purely probabilistic medium-induced modification. 

In VNI/BMS a set of partons is produced by time- like 
branching after an inelastic scattering. The branching 
process evolves the initial time-like virtuality Q of the 
outgoing parton pair downwards to a cutoff Qq ~ 1 GeV, 
converting virtuality into energy and transverse momen- 
tum of radiated partons. The formation time of these 
evolving partons Q is computed for each intermediate 
stage of the branching process. The hardest gluon in the 
final state of the time-like branching is selected to be the 
probe parton. 

The formation times 
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for all sub-branchings that occur during the time-like 
branching procedure leading up to the production of the 
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probe gluon are summed, giving the initial formation 
time for the probe gluon. The selection of the hardest 
gluon to represent the shower reflects the dominance of 
gluon re-scattering in QCD interference processes. 

The probe parton is allowed to propagate through the 
medium and interact elastically with other partons in the 
medium. After the n-th scattering its formation time is 
recalculated as 

r? = " 2- (7) 

(Xl + E i= i qj-,i) 

This tends to decrease the formation time relative to the 
original value with each additional scattering, leading to 
an increased energy loss rate. This simulates the emission 
of the shower with coherent interference from n centers. 
After the formation time expires, the radiation is consid- 
ered to have decohered from the initiating particles and 
all involved partons may interact and radiate. 

The performance of this prescription has been verified 
by considering the energy loss of a 100 GeV jet inserted 
into a medium at fixed temperature T — 0.35 GeV, see 
Fig. [2] The elastic energy loss in the VNI/BMS parton 
cascade has been previously shown to agree well with 
perturbative calculations The measured radiative 

energy loss is well described by BDMPS-Z [18] which 
gives, 

-AE= a ^ f fLHo g ±. (8) 
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By fitting the radiation only curve (the blue open 
squares in Fig. [2]) with this functional form we match the 
leading coefficient to within 20%. 

The elastic energy loss process is dominant for short 
times since the average formation time arising from the 
first hard collision of the leading-parton with the medium 
is (t/) ~ 5 fm /c. This leads to the initial suppression 
of the radiative process as we have required that no en- 
ergy be lost from the radiating parton until the formation 
times of the radiated quanta have expired. 

IV. RESULTS 

We generated candidate dijet events from Monte-Carlo 
pp collisions at y/s = 200 GeV using PYTHIA 8. Anti- 
kt reconstruction was used to identify satisfactory di- 
jet events 20 . The following kinematic requirements 
were always imposed: Pt, m in — 5 GeV, A|r7 12 | < 1.1, 
A0i2 > 7r/2 where A7712 and A0i2 are measured be- 
tween the two reconstructed jets in the event. The con- 
stituent partons of each identified jet were translated 
from PYTHIA into the parton cascade medium for mod- 
ification. The two jets in each event are evolved sep- 
arately. Table [I] shows the ranges of each parameter we 
included in the study. The study contains a total of 16000 
dijet events for each combination of parameters. For a 
given medium radius R me d we generate dijet path lengths 
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FIG. 2: The total energy loss (red closed circles) for the 
leading parton in a 100 GeV jet in a medium at 
T = 0.35 MeV, the contributions from elastic and 
radiative processes (black open circles and blue open 
boxes) are shown. A fit to ^ is shown as the black 
solid line. 



parameter 


range of values 


jet cone angle R 


0.2, 0.3, 0.4 


lead energy cut E t ,t 


25, 35, 50 GeV 


interaction mechanism 


elastic, elastic & radiation 


Temperature T 


250, 350, 450 MeV 


strong coupling a s 


0.25, 0.3, 0.35, 0.4, 0.6 


medium radius Rmed 


1, 3, 5 fm 



TABLE I: 16000 dijet events were generated for each 
combination of the above factors 



by uniformly sampling 2d jet vertices within a circle of 
radius R m ed and then generating a chord with a uniform 
angular distribution. The two chord segments define the 
path lengths for the dijet pair. The leading jet is defined 
as the jet with the highest energy at the point of obser- 
vation, the leading jet from the vacuum event may not 
be the leading jet after medium evolution. 

We consider the dijet asymmetry A,, the partonic 
jet radial profile and the partonic distribution of energy 
within the jet in terms of the variables z — pt/E t cos(Ar) 
and j t — ptsin(Ar), where p t is the transverse momen- 
tum of a parton within the jet, E t is the total transverse 
momentum of the jet and Ar = \J (A</> 2 + Ar/ 2 ) is the 
angle between this parton and the jet axis. 

We have previously applied VNI/BMS [5] with a static 
medium to jets at LHC scales y/s = 2.76 ATeV, see 
Fig. [3] Here we followed the dijet kinematics as de- 
scribed by CMS 0, E u > 120 GeV, E t , s > 50 GeV 
and A</> 12 > 27r/3. A Gaussian smearing was applied 
to reconstructed jet energies to approximate detector re- 
sponse. In this case a medium with T — 350 MeV and 
Rmed = 3 fm gives a relatively good fit to the CMS results 
PJ. This suggests that VNI/BMS can produce reasonable 
jet modification. Given the complexity of the actual de- 
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tector response we present all further results without any 
attempt at smearing with the hope that particular detec- 
tor responses can be folded into the data as required. 
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FIG. 3: The VNI/BMS post diction of the dijet 
asymmetry as measured at CMS, jets are reconstructed 
with Anti-Kt R = 0.5 and a smearing 
N(25f, \ .2\fW t ) is applied to simulate the detector 
response. 



We now consider the modification of the RHIC dijet 
asymmetry under variation of all the factors in our de- 
sign. In Fig. [4] we show the variation of Aj with the 
medium radius for two different medium temperatures. 
The higher temperature has a q roughly double the lower, 
see Fig. [I] The vacuum jet distribution falls off very 
quickly, jets with an initial E t / > 65 GeV are exceedingly 
rare. Increasing the medium temperature from 250 MeV 
to 350 MeV leads to a depletion of jets with a small mod- 
ification shifting the dijet distribution towards higher Aj 
values. By applying an energy cut to differentiate the 
leading and sub-leading jets in a dijet pair we have bi- 
ased our results. We select for leading jets which travel 
a short distance, and are less modified, along with sub- 
leading jets that travel a very long distance and so experi- 
ence a larger integrated q. This surface bias gives a large 
contribution to Aj relative to the vacuum result. For the 
remainder of the analysis we have fixed R m ed = 5 fm . 

The joint leading and sub-leading energy loss distribu- 
tion for elastic only interactions and the full simulation 
is shown in Fig. [5] for jets with T me ^ — 350 MeV and 
E t ,i > 20 GeV. In both cases the distribution is peaked 
at zero, most jets simply don't lose any energy which 
contributes to the relatively large peaks at Aj ~ 0. The 
leading partons may lose energy by interacting with the 
medium this energy has to be transported outside of the 
jet cone for the E t of the jet to be modified. The line 
AE t .e — AE t , s is also responsible for the peak at small 
Aj. The leading jet rarely loses much energy in either 
scheme while the sub-leading jet is noticeably more mod- 
ified in the radiative scheme. In Fig. [6] the distribution 
of energy loss against distance traveled is shown. The 
leading jet shows a very strong surface bias. The leading 
jets travel a few fm and lose no energy. The sub-leading 



jet distribution is also peaked at z = but this peak is 
smaller and the bulk of the distribution is spread into a 
wide range of path lengths and energy losses. An under- 
standing of the path length dependence of the leading 
and sub-leading jet energy loss is vital for understanding 
the tomographic implications of Aj. 

Since q cx T 3 the dijet asymmetry depends strongly 
upon the medium temperature. In Fig. [7] we show re- 
sults of varying the medium temperature. Increasing 
medium temperature and therefore q leads to increased 
jet-medium interactions and a strong swing in the ob- 
served Aj . The difference between the elastic only results 
and the full simulation including radiation is remarkably 
small at lower temperatures. The modification of the ra- 
diative jets is somewhat greater at T — 0.45 GeV. For 
jets at these scales the asymmetry is relatively insensitive 
to the details of the jet interactions. 

In Fig.|8]we show Aj as a function of the anti-kt cone 



angle R for jets in mediums with T = 250 MeV and 
T = 350 MeV. As R is increased the amount of di- 
jet modification is significantly reduced. It is important 
to note that the only partons which can be included as 
part of the measured jets were either directly created by 
the jet or are those which have scattered with jet par- 
tons. The thermal medium is currently artificially ex- 
cluded from the jet finder, this removes uncertainty as- 
sociated with background removal. As R increases more 
of the relatively soft radiated partons and forward scat- 
tered medium partons are included in the jet definition 
along with original hard core. This leads to higher re- 
constructed jet energies at larger i?'s which in turn leads 
to the observed reduced Aj . The effect is proportional to 
the distance traveled by the sub-leading jet, jets which 
have traveled shorter distances have built up less of a 
cloud of soft partons and picked up fewer medium par- 
tons by elastic forward scattering. The rate of transverse 
diffusion of these soft partons is proportional to q, so the 
jet-cone/cloud will be wider at higher medium tempera- 
tures. 

The leading jet energy cut is varied in Fig.[9j at fixed 
T = 350 MeV and R = 0.2, the dijet asymmetry ap- 
pears to be relatively insensitive to this parameter. We 
explore the variation of the strong coupling constant a s 



in Fig. 10 The strength of elastic interactions are scales 
with the strong coupling and so g cx c? s . The radia- 
tive process itself also depends upon the strong coupling, 
the total cross section for emission in the DLA (double 
log approximation) scheme is fijv oc {gllog 2 ) N where the 
large logs arise from the kinematics of small angle radi- 
ation [3T]. Given these considerations it is interesting to 
note that the relative similarity of the results with and 
without radiation. 

Let us now examine the modification of the radial jet 
profile under the same set of factors. We define the radial 
jet profile as the ratio of jet energy reconstructed within 
a certain jet-cone radius R relative to the reconstructed 
energy at R = 1. This gives a normalized radial profile 
which is very sensitive to variations in the medium tern- 
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FIG. 4: The dijet asymmetry as a function of the path length of the leading and sub-leading jets, for a medium with 
T = 250 MeV (left) and T = 350 MeV (right). Both figures show anti-kt reconstructed jets of radius R = 0.2 with 
Et,e > 20 GeV and all events include both radiative and elastic energy loss. 




FIG. 5: The joint energy loss distributions for leading and sub-leading jets at T me d — 350 MeV with 
: 5 fmand anti-kt R = 0.2. The left figure shows the elastic only mode, right figure shows elastic + radiation. 
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FIG. 6: Showing the distribution of energy loss against path length for leading (left) and sub-leading (right) jets at 
Tmed = 350 MeV with R m ed = 5 fm and anti-kt R = 0.2 in the full radiation+elastic simulation. 
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FIG. 7: Showing the influence of the medium temperature and the jet-medium interaction mechanism for dijets with 
Eti > 20 GeV. The left panel shows elastic scattering only, the right panel shows the full elastic and radiative 
process. All events are for anti-kt reconstructed jets with R = 0.2 and R m ed — 5 fm . 
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FIG. 8: Showing the variation of Aj under variation of R the Anti-Kt jet cone radius, elastic only (left) and elastic 
+ rad (right). The medium radius is kept fixed at R me d = 5 fmwith T = 0.35 GeV and E t ^ > 20 GeV. 
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perature, the jet interaction mechanism and the strong 
coupling. This observable could be experimentally ob- 
tained by a similar process of iterative jet reconstruction. 
In Fig.[TT]we show the radial profile at two medium tem- 
peratures for both elastic and radiatively modified jets. 
The leading jet profiles (solid lines) are somewhat mod- 
ified compared to the vacuum jets (black profiles), we 
see very little separation between the radiative and elas- 
tic only leading jet profiles. The sub-leading jet profiles 
are dramatically modified compared to the vacuum and 
leading jet profiles. At both temperatures the elastic and 
radiative sub-leading jet profiles clearly separate, the ra- 
diative sub-leading jets become broader and softer than 
the elastic only. Both sets of sub-leading jets become 
much broader and softer compared to the leading jets. 
This is in line with our observations of the energy loss 
distributions as shown in Fig.[5j In Fig. 12 we show the 
jet profiles' sensitivity to the strong coupling. The lead- 
ing jet profiles are somewhat modified at a s — 0.6, the 
sub-leading jet profiles soften smoothly with increasing 
a s . It is interesting to note that the shape of these soft 
sub-leading profiles is different for elastic (left) and radia- 
tive (right) jets, the elastic profiles become increasingly 
convex as a s — > 0.6 while the radiative profiles remain 
concave. 

Finally we examine the partonic fragmentation distri- 
butions in the longitudinal and transverse variables jt 
and z. Examining the single inclusive jet distribution at 
RHIC scales yfs = 200 GeV and E t ~ 15 - 65 GeV we 
find that the average partonic multiplicity in a vacuum 
PYTHIA jet is (N parton ) = 5.6 which is roughly tripled 
during hadronization (N nac [ ron } — 14.6, this hadroniza- 
tion process will conserve the total E t of the the jet but 
in so doing will redistribute the components of the mo- 
mentum transverse to the jet axis. This means that the 
partonic and j t and z distributions cannot be directly 
compared with experiment. What may be a surprise is 
the observed large peak in the z distribution at 1 cor- 
responding to jets with a single parton, this is seen in 
the vacuum and modified jet results. This peak is a 
consequence of the vacuum radiation evolution which is 
kinematically limited at RHIC regimes. We note that at 
LHC jet scales E t ~ 50 - 300 GeV and = 2.76 TeV 
{Nparton) = 19.5 and (Nhadron) = 46.5, at these scales 
the z = 1 peak vanishes. 



In Fig.[13] we show the VNI/BMS partonic z = 
Pt/Ef cos(Ar) distributions for elastic and radiative jets 
alongside the vacuum results. Note the strong peak at 
z = 1 in both cases, this peak is actually enhanced by 
the evolution. In the elastic case this enhancement comes 
from a depletion of the low z distribution. If a jet ini- 
tially contains only one or two soft partons and a hard 
core the soft partons will most likely be scattered out of 
the jet cone. This leads to an enhancement of jets mea- 
sured with z = 1 and a depletion of the low end of the 
spectrum. The same process takes place in the radiative 
case, the right panel of Fig. [13] but here the soft radiated 
partons tend to fill the small z distribution while the mid 



z range sees the most depletion. These middle z partons 
interact with the medium producing a shower of soft ra- 
diation. The enhancement of the z — 1 peak is smaller in 
the radiative case than in the elastic case as these hard 
core partons radiate and soften the jets. As the jet cone 
radius is increased the peak at z = 1 is diminished as 
more of the soft scattered partons are recaptured by the 
jet definition this also enhances the peak at small z in the 
radiative case, see Fig. 14 The transverse jt = Pt sin(Ar) 



The transverse jt 
distribution is shown in Fig.[l5] for elastic and radiative 
jets. The evolved j t distribution is softened in both cases 
as interactions with the medium and radiation add more 
soft partons to the jet. The radiative jets have a slightly 
steeper profile than the elastic only jets, radiated soft 
partons are more strongly confined to small j t . 



V. CONCLUSIONS 

We have examined the modification of RHIC scale di- 
jets by a pQCD medium over a range of temperatures 
and coupling constants. We find that the dijet asymme- 
try is quite sensitive to the amount of medium encoun- 
tered by the jets, and very sensitive to variations in q. We 
do not observe a strong dependence of Aj upon the jet 
interaction mechanism which is intriguing. The leading 
jets show a strong surface bias and lose very little energy 
making their E t a reasonable proxy for their initial en- 
ergy. As such Aj appears to give a reasonable measure 
of the modification of sub-leading jets at these scales. 

The jet profiles we present show the strong modifi- 
cation of the sub-leading jets as a function of temper- 
ature and strong coupling. They also differentiate be- 
tween interaction mechanisms, understanding the nature 
of the jet-medium interaction as well as its strength is of 
paramount import. Although these profiles may well be 
sensitive to the non-trivial hadronization process, we be- 
lieve that they could provide a very informative window 
into the real nature of jet interactions with the QGP. 

The fragmentation distributions provide another in- 
teresting window into jet modification. The observed 
structures in z reveal the importance of the interplay 
of vacuum radiation and the fragmentation process. Al- 
though hadronization will certainly change these distri- 
butions these results already demonstrate the importance 
of medium induced transverse diffusion in the jet modi- 
fication process. 

Dijets at RHIC scales are likely to be strongly modified 
by the presence of the deconfined QGP medium. The ob- 
servables we have discussed are sensitive to many aspects 
of this modification and suggest that further jet measure- 
ments at RHIC will provide valuable insights into the na- 
ture of the QGP and into the applicability of pQCD jet 
suppression models . 
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FIG. 9: The influence of applying cuts to the leading jet energy on A, . The left panel shows only elastic scattering 
processes, the right panel shows radiation and elastic scattering. All events shown are for anti-kt reconstructed jets 

with R = 0.2 and R m ed — 5 fm . 
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T = 250 MeV reconstructed with anti-kt R = 0.2. 
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FIG. 11: Variation of the jet radial profile at T = 250 MeV (left) and T = 350 MeV (right). Leading jet profiles are 
shown as solid lines and sub-leading jets are dashed, elastic only jets are red and radiative + elastic jets are orange. 
Note the strong softening of the sub-leading jet profiles and the clear separation of the elastic and radiative cases. 
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FIG. 12: Variation of the jet radial profile at T = 350 MeV. Leading jet profiles are shown as solid lines and 
sub-leading jets are dashed, elastic only jets are shown on the left and radiative + clastic jets are on the right. 
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FIG. 13: The z distribution for partonic jets at T = 350 MeV for jets reconstructed at R — 0.2 for elastic jets (left) 
and radiative+elastic (right). Leading jets are shown as open triangles and subdeading jets are shown as crosses. 
The vacuum distributions are shown as red solid lines and the evolved distributions as blue dashed lines. 
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and radiative+elastic jets (right). Leading jets are shown as open triangles and subdeading jets are shown as 
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